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ABSTRACT 

We study the generation of large-scale vortices in rotating turbulent convection by means of Cartesian direct 
numerical simulations. We find that for sufficiently rapid rotation, cyclonic structures on a scale large in 
comparison to that of the convective eddies, emerge, provided that the fluid Reynolds number exceeds a critical 
value. For slower rotation, cold cyclonic vortices are preferred, whereas for rapid rotation, warm anti-cyclonic 
vortices are favoured. In some runs in the intermediate regime both types of cyclones co-exist for thousands 
of convective turnover times. The temperature contrast between the vortices and the surrounding atmosphere 
is of the order of five per cent. We relate the simulation results to observations of rapidly rotating late-type 
stars that are known to exhibit large high-latitude spots from Doppler imaging. In many cases, cool spots are 
accompanied with spotted regions with temperatures higher than the average. In this paper, we investigate a 
scenario according to which the spots observed in the temperature maps could have a non-magnetic origin due 
to large-scale vortices in the convection zones of the stars. 
Subject headings: Hydrodynamics - convection - turbulence 



1. INTRODUCTION 

Rotating turbulent convection is considered to play a crucial 
role in the generation of large-scale magnetic fields (Moffatt 
1978; Krause & Radler 1980; Rudiger & Hollerbach 2004) 
and differential rotation of stars (Rudiger 1989). The interac- 
tion of rotation and inhomogeneous turbulence leads to the so- 
called a-effect, which can sustain large-scale magnetic fields 
(e.g. Brandenburg 2001; Kapyla et al. 2009). However, in 
many astrophysically relevant cases large-scale shear flows 
are also present, which further facilitate dynamo action by 
lowering the relevant critical dynamo number. In the Sun, for 
example, the entire convection zone is rotating differentially 
(cf. Schou et al. 1998; Thompson et al. 2003), and a merid- 
ional flow towards the poles is observed in the near surface 
layers (e.g. Zhao & Kosovichev 2004). These flows are most 
often attributed to rotationally influenced turbulent angular 
momentum and heat transport (cf. Rudiger 1989; Robinson 
& Chan 2001; Miesch et al. 2006; Kapyla et al. 201 lb). In the 
solar case the large-scale flows and also the magnetic activity 
are largely axisymmetric (e.g. Pelt et al. 2006). This means 
that the sunspots, which are concentrations of strong magnetic 
fields, are almost uniformly distributed in longitude over the 
solar surface. The fact that we observe the sunspots and can 
attribute magnetic fields to them, has strongly influenced the 
interpretation of data from stars other than the Sun. 

The giant planets Jupiter and Saturn are also likely to have 
outer convection zones (e.g. Busse 1976), but they rotate 
much faster than the Sun. Bands of slower and faster rotation 
alternate in their atmospheres, reminiscent of rapidly rotating 
convection (e.g. Busse 1994; Heimpel & Aurnou 2007). How- 
ever, especially in Jupiter, large spots in the form of immense 
storms are observed (Marcus 1993). Remarkably, the largest 
of these, the Great Red Spot, has persisted at least 180 years. 
Similar features are observed also in Saturn (e.g. Sanchez- 
Lavega et al. 1991) and other giant planets. The spots on gi- 
ant planets are not of magnetic origin although dynamos are 
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likely to be present in the interiors of the planets. Thus their 
explanation is probably related to hydrodynamical processes 
within the convectively unstable layers. 

Late-type stars with higher rotation velocities in compari- 
son to the Sun, on the other hand, often exhibit light curve 
variations that are usually interpreted as large spots on the 
stellar surface (e.g. Chugainov 1966; Henry et al. 1995). In 
some cases the observational data can be fitted with a model 
with two large spots at a 180 degree separation in longitude 
Berdyugina & Tuominen (1998). There is also evidence that 
these 'active longitudes' are not equal in strength (e.g. Lehti- 
nen et al. 2011; Lindborg et al. 2011), and that the relative 
strenght of the spots can, at least temporarily, reverse in a pro- 
cess dubbed 'flip-flop' (cf. Jetsu et al. 1993). One interpreta- 
tion of the data is that the spots are of magnetic origin and that 
the flip-flops are related to magnetic cycles reminescent of the 
solar cycle (e.g. Berdyugina et al. 1998). On the other hand, 
it has been proposed that the flip-flops are only short-term 
changes related to the activity cycle, while the structure gen- 
erating the temperature minima would migrate in the orbital 
reference frame, that could be interpreted as an azimuthal dy- 
namo wave (e.g. Lehtinen et al. 2011; Lindborg et al. 2011). 
Again, this interpretation relies on the magnetic nature of the 
cool spots. 

The cool spots detected by photometry and Doppler imag- 
ing using spectroscopic observations have been taken as an in- 
direct proxy of the magnetic field on the stellar surface, deriv- 
ing from the analogy to sunspots - strong magnetic field hin- 
ders convection and causes the magnetized region to be cooler 
than its surrounding. Zeeman-Doppler imaging of spectropo- 
larimetric observations (e.g. Semel 1989; Donati et al. 1989; 
Piskunov & Kochukhov 2002; Carroll 2007) provides means 
to directly measure the magnetic field strength and orientation 
on the stellar surface. In the study of Donati et al. (1997) spec- 
tropolarimetric observations of several stars were collected 
during 23 nights extending over a five year interval. They 
report that the Zeeman signitures of the cool stars almost al- 
ways exhibit a very complex shape with many successive sign 
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reversals. This points to a rather complicated field structure 
with different magnetic regions of opposite polarities. Fur- 
thermore, the magnetic regions detected were mostly 500 to 
1,000 K cooler than, and sometimes at the same temperature 
as, but never warmer than the surrounding photosphere. In the 
published temperature and magnetic field maps for AB Dor 
(Donati & Collier Cameron 1997), however, no clear corre- 
lation between temperature and magnetic field strength can 
be seen: in the temperature maps a pronounced cool polar cap 
with weak fringes towards lower latitudes are visible, whereas 
the strongest magnetic fields are seen as patchy structures at 
lower latitudes with a clearly different distribution than the 
temperature structures. Similar decorrelation of temperature 
minima and magnetic field strength has been reported with 
the same method for different objects (e.g. Donati 1999; Jef- 
fers et al. 2011), and also for the same objects with different 
methods (e.g. Hussain et al. 2000; Kochukhov et al. 2011). 
The phenomenon, therefore, seems to be wide-spread, and 
method-independent. 

One possible explanation to the decorrelation of magnetic 
field and temperature structures could be that there is sim- 
ply less light coming from the spotted parts than from the 
unspotted surface. Thus the Zeeman signatures from cool 
spots may be "drowned" in the signal from the unspotted sur- 
face or bright features. However, this should lead to system- 
atic effects where the detected magnetic field strength would 
be correlated with the surface temperature. The least sqau- 
res deconvolution technique (LSD, e.g. Donati et al. 1997), 
which is necessary for enhancing the Zeeman signal, may in- 
fluence the temperature and magnetic Doppler imaging differ- 
ently. The latitudes of any surface features in Doppler images 
are always more unrealiable than the longitudes, a fact that 
will not make a comparison of temperature and magnetic field 
maps any easier. One could thus expect, that there could be 
artificial discrepancies in the latitudes of magnetic and tem- 
perature features. Still, the lack of connection between even 
the longitudes of cool spots and magnetic features is surpris- 
ing. 

In this paper we consider a completely different scenario, 
according to which the formation of temperature anomalies 
on the surfaces of rapidly rotating late-type stars could oc- 
cur due to a hydrodynamical instability creating large-scale 
vortices, analogously to the giant planets in the solar system. 
To manifest this mechanism in action, we simulate rotating 
turbulent convection in local Cartesian domains, representing 
parts of the stratified steUar convection zones located near the 
polar regions. We show that under such a setting, large-scale 
vortices or cyclones are indeed generated provided that the ro- 
tation is sufficiently rapid and the Reynolds number exceeds 
a critical value. Depending on the handedness of the vortex, 
which on the other hand depends on the rotation rate, the re- 
sulting spot can be cooler or warmer than the surrounding at- 
mosphere. 

We acknowledge that our model is rather primitive, lacking 
realistic radiation transport, spherical geometry, and relying 
on a poly tropic setup for the stratification so detailed compar- 
ison with observations is not possible at this point. However, 
the main purpose of the present paper is to show a proof of 
concept of the existence of large-scale vortices with tempera- 
ture anomalies close to those observed in rapidly rotating hy- 
drodynamic convection. We also note that similar large-scale 
cyclonic structures have recently been reported from large- 
eddy simulations of turbulent convection (Chan 2003, 2007). 
We make comparisons to these studies when possible. 



2. THE MODEL 

Our model setup is similar to that used by Kapyla et al. 
(2009) but without magnetic fields. A rectangular portion of 
a star is modeled by a box situated at colatitude 6. The box 
is divided into three layers: an upper cooling layer, a convec- 
tively unstable layer, and a stable overshoot layer (see below). 
We solve the following set of equations for compressible hy- 
drodynamics: 
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where D/Dt = d/dt + U • V is the advective time derivative, 
y is the kinematic viscosity, K is the heat conductivity, p is 
the density, U is the velocity, g = -gz is the gravitational 
acceleration, and £1 - Qo(- sin 6, 0, cos 6) is the rotation vec- 
tor. The fluid obeys an ideal gas law p = (y - l)pe, where 
p and e are pressure and internal energy, respectively, and 
7 = cp/cy = 5/3 is the ratio of specific heats at constant 
pressure and volume, respectively. The specific internal en- 
ergy per unit mass is related to the temperature via e - cyT. 
The rate of strain tensor S is given by 



Sij=^iU,j + Uj,d- 



i<5,-,V • U. 



(4) 



The last term of Eq. (3) describes cooling at the top of the do- 
main. Here t(z) is a cooUng time which has a profile smoothly 
connecting the upper cooling layer and the convectively un- 
stable layer below, where t — > oo. 

The positions of the bottom of the box, bottom and top of 
the convectively unstable layer, and the top of the box, re- 
spectively, are given by (21,22,23,24) = (-0.85,0, 1, l.l5)d, 
where d is the depth of the convectively unstable layer. Ini- 
tially the stratification is piecewise polytropic with polytropic 
indices (mi,m2,m3) = (3, 1, 1), which leads to a convectively 
unstable layer above a stable layer at the bottom of the do- 
main. In a system set up this way, convection transports 20 
per cent of the total flux (cf. Brandenburg et al. 2005). Due 
to the presence of the cooling term, a stably stratified isother- 
mal layer is formed at the top. The horizontal extent of the 
box, Lii = Lx - Ly, is 4d. All simulations with rotation are 
made at the North pole, corresponding to = 0°. The simu- 
lations were performed with the Pencil Code\ which is a 
high-order finite difference method for solving the compress- 
ible equations of magnetohydrodynamics. 

2. 1 . Units and nondimensional parameters 
Nondimensional quantities are obtained by setting 

^ = 5 = Po = cp = 1 , (5) 

where po is the initial density at 22- The units of length, time, 
velocity, density, and entropy are 

[x\ = d, m= yfdil, m = yfd^, 

[p] = Po , [s] = cp. (6) 
' http://code.google.coin/p/pencil-code/ 
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We define the Prandtl number and the Rayleigh number as 

where xo - K/ (pmCp) is the thermal diffusivity, and Pm is the 
density in the middle of the unstable layer, Zm - 3(^3 - zi)- 



Pr=^ 



1 d,s\ 
CP dz/o 



The entropy gradient, measured at Zn 
hydrostatic state, is given by 



in the nonconvecting 



1 ds\ 
CP dzJo 



V-Va. 

Hp 
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where V - V^tj is the superadiabatic temperature gradient with 
Vad - 1 - 1 /y, V = {dlnT/dln p)r^ , and where Hp is the pres- 
sure scale height. The amount of stratification is determined 
by the parameter ^0 - iy - ^)e\il{gd), which is the pressure 
scale height at the top of the domain normalized by the depth 
of the unstable layer. We use ^0 - 1 /3 in all cases, which re- 
sults in a density contrast of about 23 across the domain. We 
define the Reynolds and Peclet numbers via 



Re 



Pe 



Pr Re 



(9) 



vkf xokf 

where kf - In/d is adopted as an estimate for the wavenum- 

ber of the energy-carrying eddies, and Mrms - yj^u^. This 

definition neglects the contributions from the large-scale vor- 
tices that are generated in the rapid rotation regime. Note that 
with our definitions Re and Pe are smaller than the usual one 
by a factor In. The amount of rotation is quantified by the 
CorioUs number, defined as 



Co = 



2Qo 



(10) 



We also quote the value of the Taylor number, 

Ta = (2Qo^/'/vf , (11) 
which is related to the Ekman number via Ek = Ta"'^^. 

2.2. Boundary conditions 

The horizontal boundaries are periodic for all variables. 
Stress-free conditions are used for the velocity at the vertical 
boundaries. 



U,,, = Uy,, = U, = 0. 



(12) 



Temperature is kept constant on the upper boundary and the 
temperature gradient 



dT 

dz 



cv(r- l)(m+ 1) 



(13) 



is held constant at the lower boundary, yielding a constant 
heat flux Fq = -KdT/dz through the lower boundary. 

3. RESULTS 

We perform a number of numerical experiments in order to 
determine the conditions under which large-scale cyclones are 
excited. The basic input parameters and some key diagnostic 
outputs of the simulations are listed in Table 1. We perform 
a few (Set A) or a single (Sets B, C, and D) progenitor run 
with a given Peclet number in each Set from which the rest of 
the runs are obtained by continuing from a saturated snapshot 
and changing the value of the kinematic viscosity v in order to 





Figure 1. Upper panel: total velocity from Runs A5 and A7. Lower panel: 
velocity components (black), ^ji^. (red), and ^hF (blue) from Run A7. 
The jump at ftinns^^f ~ 500 is due to a lowering of v at this point. 

change Re. The higher resolution run D2 was remeshed from 
a lower resolution case Dl. 

3.1. Excitation of large-scale vortices 

We perform several sets of runs where the Peclet number 
and input energy flux are constant, whereas the Reynolds and 
Coriolis numbers are varied. We are limited to exploring a 
small number of cases due to the slow growth of the vortices, 
see Table 1 . Typically the time needed for the saturation of the 
cyclones is several thousand convective turnover times (see 
Fig. 1). Thus many of our runs were ran until the presence or 
the absence of the cyclones was apparent. 

We find that a reliable diagnostic indicating the presence 
of large-scale vortices is to compare the rms-value of the to- 
tal velocity, t/rms. and the volume average of the quantity 

Mims = The latter neglects the horizontal velocity 

components which significantly grow when large-scale cy- 
clones are present (see the lower panel of Fig. 1). In the 
cyclone-free regime, irrespective of the rotation rate, we find 
that f/rms ~ Mrins Suggesting that the flow is only weakly 
anisotropic (see Table 1). In the growth phase of the vortices 
one of the horizontal velocity components is always stronger, 
but the relative strength of the components changes as a func- 
tion of time (see the lower panel of Fig. 1). This undulation is 
related to quasi-periodic changes of the large-scale pattern of 
the flow, although their ultimate cause is not clear. 

Another quantitative diagnostic is to monitor the power 
spectrum of the flow from a horizontal plane within the con- 
vection zone. A typical example is shown in Fig. 2 where 
power spectra of the velocity from the middle of the convec- 
tion zone at two different times from Run B3 are shown. The 
snapshot from tUi-m^kf = 1830 is the initial state for Run B3, 
taken from Run Bl, showing no cyclones. The power spec- 
trum shows a maximum at k/ki - 7, indicating that most of 
the energy is contained in structures having a size typical of 
the convective eddies. However, as the run is continued fur- 
ther, a large-scale contribution due to the appearance of the 
vortices, peaking at k/ki - 1 grows, and ultimately domi- 
nates the power spectrum. We note that this run was not ran 
until saturation so the peak at k/ki - 1 is likely to be even 
higher in the final state. The presence of the vortices is also 
clear by visual inspection of the flow. A typical example is 
shown in Fig. 3, where the vertical velocity component, u^, is 
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Table 1 

Summary of the runs. Here, Ma = Urms/igd)^^' and Ma^ = Jimis/C^'^)''^- Brackets indicate that the simulation has not been run upto saturation. The 
dimensionless input heat flux at the lower boundary of the box is given by Fq = Fo/(pcl). The last column indicates the presence of cyclonic (C), anti-cyclonic 

(A), or both types (A+C) of vortices. 



Run 


grid 


Ma 


Ma, 


Re 


Pe 


Pr 


Ra 


Co 


Ta 


Fo 


Cyclones 


Al 


256^ X 128 


0.048 


0.020 


33 


8 


0.24 


2.0 


lO*- 


15.5 


4.0 ■ 10** 


1.7 


10-5 


yes (A) 


A2 


256^ X 128 


0.018 


0.017 


13 


6 


0.48 


1.0 


lO'' 


14.4 


5.6 ■ lO' 


1.7 


10-5 


no 


A3 


256 X 128 


0.022 


0.019 


21 


7 


0.36 


1.3 


10*' 


12.3 


1.0 ■ 10** 


1.7 


10-5 


no 


A4 


256^ X 128 


(0.063) 


0.023 


37 


9 


0.24 


2.0 


lO*" 


10.3 


2.3 ■ 10** 


1.7 


10-5 


yes (A) 


A5 


256^ X 128 


0.021 


0.020 


16 


9 


0.48 


1.0 


lO*- 


7.9 


2.5 ■ lO' 


1.7 


10-5 


no 


AO 


o c^2 w 1 oo 
ZjO X Izo 


0.024 


0.023 


24 


9 


0.36 


1.3 


10*' 


7.0 


4.4 ■ lU 


1.7 


10-5 


no 


A7 


256 X 128 


(0.093) 


0.026 


42 


10 


0.24 


2.0 


10*' 


6.1 


1.0 ■ 10* 


1.7 


10-5 


yes (A+C) 




ZjO X IZo 


U.UZo 


A 097 
U.UZ / 


zo 


1 1 
1 1 


yi.jo 


1.3 


lO* 


J.O 


1 .o ■ lU 


1.7 


10-5 


no 


A9 


256^ X 128 


0.082 


0.028 


45 


11 


0.24 


2.0 


10*^ 


3.4 


3.6 ■ 10' 


1.7 


10-5 


yes (C) 


A9b 


256^ X 128 


(0.070) 


(0.031) 


49 


12 


0.24 


2.0 


10** 


2.1 


1.6- 10' 


1.7 


10-5 


decay 


AlO 


256^ X 128 


0.032 


0.033 


53 


13 


0.24 


2.0 


10* 


1.0 


4.0 ■ 10*^ 


1.7 


10-5 


no 


All 


256^ X 128 


0.038 


0.038 


61 


15 


0.24 


2.0 


10* 








1.7 


10-5 


no 


Bl 


256^ X 128 


0.017 


0.016 


26 


13 


0.48 


4.0 


10" 


9.7 


1.0- 10** 


8.6 


10-" 


no 


B2 


256^ X 128 


(0.021) 


(0.017) 


37 


13 


0.36 


5.4 


10'' 


9.1 


1.8 ■ 10** 


8.6 


lo-*" 


yes (A+C) 


B3 


256^ X 128 


(0.034) 


(0.020) 


63 


15 


0.24 


8.0 


lO*- 


8.0 


4.0 ■ 10** 


8.6 


K)-'' 


yes (A+C) 


CI 


256^ X 128 


0.011 


0.011 


17 


16 


0.96 


8.0 


10'' 


14.8 


1.0- 10** 


4.3 




no 


C2 


256^ X 128 


(0.014) 


(0.012) 


25 


18 


0.72 


1.1 


10^ 


13.6 


1.8 ■ 10** 


4.3 


K)-'' 


no 


C3 


256^ X 128 


(0.022) 


(0.014) 


44 


21 


0.48 


1.6 


10^ 


11.6 


4.0 ■ 10** 


4.3 


K)-'' 


yes (A) 


Dl 


256^ X 128 


0.013 


0.013 


42 


51 


1.20 


4.0 


10' 


7.2 


1.4- 10** 


1.7 




no 


D2 


512^ x256 


(0.038) 


(0.013) 


101 


49 


0.48 


1.0 


10^ 


7.5 


9.0 ■ 10** 


1.7 




yes (A+C) 




Run D3 

(=26617,0 



10 



100 



Figure 2. Power spectra of velocity from early (dashed line) and late (solid) 
times from Run B3. 



shown from the periphery of the domain for Run D2. 

The data in Table 1 suggests that large-scale vorticity is ex- 
cited provided that the Reynolds number exceeds a critical 
value, Rec. For Pe 10 (Set A) we find that Rec is around 30, 
although the sparse coverage of the parameter range does not 
allow a very precise estimate to be made. We find a similar 
value for Rcc in Sets B and C, whereas for Pe w 50 in Set D, 
the critical Reynolds number is greater than 42. In Set C, 
Runs CI and C2 were started from a snapshot of Run C3 at a 
time when vortices were already clearly developing. In both 
cases we find that the cyclones decay, suggesting that their 
presence is not strongly dependent on the history of the run. 

The critical Coriolis number in Set A is somewhere be- 
tween 2.1 and 3.4. Again a very precise determination cannot 
be made, but continuing from a saturated snapshot of Run A9 
with a somewhat lower rotation rate indicates that the vortices 
decay (Run A9b). We have limited the present study to the 
North pole {6 - 0), but vortices are also excited at least down 
to latitude = 45° in the study of Chan (2007). 




Figure 3. Vertical velocity component U- at the periphery of the box from 
Run D2. See also http://www.helsinki.fi/~kapyla/movies.html. The top and 
bottom panels show slices near the top and bottom of the convectively unsta- 
ble layer, respectively. 

3.2. Thermal properties of the cyclones 

In order to study the possible observable and other effects 
of the vortices, we ran a few simulations in Set A (Runs Al, 
A4, A7, A9, AlO, and All) to full saturation. Figures 4 and 
5 show the vertical velocity and temperature in the saturated 
regime from the six runs listed above. In the non- and slowly 
rotating cases (the two rightmost panels in the lower rows of 
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Figure 4. From left to right: vertical velocity from Zm for Runs Al, A4, A7 (upper row), A9, AlO, and All (lower row). The rotational influence is decreasing 
from left to right and top to bottom. 



Figs 4 and 5), convection shows a typical cellular pattern. Vor- 
ticity is generated at small scales at the vertices of the convec- 
tion cells, but no large-scale pattern arises. We note that long- 
lived large-scale circulation can also emerge in non-rotating 
convection (e.g. Bukai et al. 2009). However, such structures 
are not likely to be of relevance in rapidly rotating stars. 

When the rotation is increased to Co » 3.3, a cyclonic, i.e. 
rotating in the same sense as the overall rotation of the star, 
vortex appears (the lower left panels of Figs. 4 and 5). Vertical 
motions are suppressed within the vortex and it appears as a 
cool spot in the temperature slice. Increasing rotation further 
to Co K 6, also an anti-cyclonic, i.e. rotating against the over- 
all fluid rotation, warm vortex appears (the rightmost upper 
panels of Figs. 4 and 5). In Run A7 the two vortices coexist for 
thousands of convective turnover times. In the most rapidly 
rotating cases Al and A4 (the two leftmost panels in the up- 
per rows of Figs. 4 and 5) a single anti-cyclonic vortex persists 
in the saturated regime. A similar behaviour as a function of 
rotation was found by Chan (2007) from large-eddy simula- 
tions. The anti-cyclonic vortices show vigorous convection 
whereas in the surrounding regions convection appears sup- 
pressed. Due to the enhanced energy transport by convection, 
the anti-cyclones appear as warmer structures than their sur- 
roundings in the temperature slices. 

Figure 6 shows that in Runs A9 and Al the flow is in 
geostrophic balance, i.e. that the flow follows the isocontours 
of pressure for both types of vortices. The cyclone in Run A9 
shows as a low pressure area, similarly to the cyclones in the 
atmosphere of the Earth, whereas the anti-cyclone in Run Al 
coincides with a high pressure region. A weaker high pressure 
region is present also in Run A9. It is not clear whether this 
kind of single or two spot configuration lasts if the domain is 
larger in the horizontal directions, or whether a greater num- 
ber of spots appear. We find that the temperature contrast be- 
tween the spot and the surrounding medium is of the order of 



five per cent (Fig. 7) for both types of vortices. Although the 
relative temperature contrast between the vortex and the sur- 
rounding vortex-free convection seems to be a robust feature 
in the simulations, we must remain cautious when compar- 
ing the results with observations. This is due to the rather 
primitive nature of the simulations that lack realistic radiation 
transport. Convection in our model is also fairly inefficient by 
design, only 20 percent of the total flux being carried by it. 

3.3. Dynamo considerations and discussion 

Figure 8 shows the horizontally averaged kinetic helicity, 
CO ■ u, where co -V y.u, from Runs A9 and Al from the ini- 
tial, purely convective cyclone-free, and final fully saturated 
stages of the simulations. The data is averaged over a period 
of roughly 60 convective turnover times in each case. We find 
that in Run A9, where a cool cyclonic vortex appears, there 
is almost no change in the kinetic helicity between the initial 
and final stages of the simulation. In this run convection, and 
thus vertical motions, are largely suppressed within the vor- 
tex (see Fig. 4). Furthermore, the dominant contribution to 
the vorticity due to the cyclone arises via the vertical compo- 
nent 0), - dxUy - dyUx, which is positive for a cyclonic vortex. 
These two effects seem to compensate each other and the he- 
licity within the cyclone is not greatly enhanced or depressed 
with respect to the surroundings. This would indicate that the 
influence of the cyclonic vortices on the magnetic field am- 
plification would be minor, as the helicity remains unaltered. 
On the other hand, the strong horizontal motions connected 
to the cyclone might be able to amphfy the field by advecting 
the field lines. 

In Run Al, on the other hand, a more pronounced effect is 
seen, and the helicity is decreased up to a factor of two in the 
saturated stage (see the right panel of Fig. 8). This change is 
brought about by the different handedness of the vorticity in 
the anti-cyclone and by the vigorous convection within it (see 
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Figure 5. From left to right: temperature from Zm for Runs Al, A4, A7 (upper row), A9, AlO, and Al 1 (lower row). The rotational influence is decreasing from 
left to right and top to bottom. The solid and horizontal lines on the leftmost panels correspond to line plots shown in Fig. 7. 

and Dl of KKB09) are clearly in the parameter regime ex- 
ceeding the critical values found here, and another four runs 
(A5, A6, B5, and CI of KKB09) where the parameters were 
close to marginal. The Reynolds and Coriolis numbers for 
these runs were calculated from the saturated state of the dy- 
namo which in all cases lowers the turbulent velocities some- 
what, decreasing the Reynolds and increasing the Coriolis 
numbers correspondingly. Furthermore, a different defini- 
tion of the Reynolds number was used by KKB09 as in the 
present study. A reanalysis of the data of KKB09 suggests 
that early stages of cyclone formation are in progress in all 
of the runs listed above. However, the magnetic field grows 
on a significantly shorter timescale than the cyclones, and the 
magnetic field saturates akeady before thousand convective 
turnover times. None of the runs was continued much fur- 
ther than twice that, making it impossible to decide in favor 
or against the maintenance of vortices based on these runs. 

Nonetheless, indications of growing cyclones appear in the 
kinematic regime, i.e. when the magnetic field is weak in 
comparison to kinetic energy of the turbulence, but they are 
far less clear, or even absent, when the magnetic field satu- 
rates. This raises two related questions: firstly, are the vortices 
responsible for the emergence of the large-scale magnetic 
fields, and secondly, can the vortices coexist in the regime 
where strong magnetic fields are present? The current data 
suggests that the presence of the vortices is not essential for 
the large-scale magnetic fields which persist throughout the 
saturated state, whereas the vortices remain less prominent 
or suppressed. This fact is related to the second issue. As 
noted above, the simulations of KKB09 are too short for the 
vortices to fully saturate. Thus, we cannot conclusively say 
whether the lack of the vortices in the dynamo regime is due 
to the magnetic field simply reducing the Reynolds number 
below the critical value, or a direct influence of the Lorentz 
force on the growing vortices. We will address the questions 



Figure 6. Pressure (colors) and horizontal flows from the middle of the con- 
vection zone in Runs A9 (left panel) and Al (right panel). 

the upper row of Fig. 4). The combination of these produces 
significantly greater helicity in the anti-cyclones, but a pre- 
dominantly different sign than the surroundings and leads to 
an overall decrease noted in Fig. 8. The decreased amount of 
helicity would indicate weaker amplification of the magnetic 
field by anti-cyclones compared to their surroundings. Again, 
the strong horizontal motions might counteract by amplifying 
the field by advection. 

The simulations presented here were performed with a 
setup identical to that used in (Kapyla et al. 2009, hereafter 
KKB09) to study large-scale dynamo (LSD) action in rotating 
convection. In KKB09 the generation of large-scale magnetic 
fields, given that the Coriolis and magnetic Reynolds numbers 
exceeded critical values, were reported. The critical Coriolis 
number for LSD action was found to be roughly four, which 
is close to the critical value for the cyclones to emerge. 

The relation of the two phenomena is an interesting ques- 
tion, that can be only partially answered by the existing mag- 
netohydrodynamic runs. This is because the fluid Reynolds 
number in the runs of KKB09 was in most cases lower than 
Rec required for the vortices to appear. Only two runs (AlO 
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Figure 7. Temperature as a function of x from a quiescent (solid lines) and 
cyclonic (dashed) regions for Runs A9 (left panel) and Al (right panel). The 
positions of the cuts are indicated in the leftmost panels of Fig. 5 with corre- 
sponding linestyles. The normalization factor T is the horizontal average of 
the temperature. 




Figure 8. Horizontally averaged kinetic helicity w ■ u as a function of z from 
a quiescent (solid lines) and cyclonic (dashed) states for Runs A9 (left panel) 
and Al (right panel). The vertical dotted lines at z = and z = d indicate the 
bottom and top of the convectively unstable layer, respectively. 

related to magnetic fields and dynamo action in more detail in 
a forthcoming paper. 

3.4. Observational implications 

If large-scale cyclones such as those found in the present 
study occur in real stars, they will cause observational signa- 
tures on the stellar surface due to their lower or higher tem- 
perature. The temperature contrasts seen in the surface maps 
derived by Doppler imaging are somewhat stronger than the 
value of roughly five percent found in this study; for instance, 
on the surface of the active RS CVn binary 11 Peg, analysed 
by Lindborg et al. (2011) and Hackman et al. (2QIY), the 
coolest spot temperatures, depending on the season, are 10- 
20 per cent below the mean surface temperature. Similar spot 
temperatures have also been obtained by analysing molecu- 
lar absorption bands, but cooler stars seem to have a lower 
spot contrast (O'Neal et al. 1998). Taken that the numeri- 
cal model is quite simple, for instance in the sense that the 
transport of energy by convection is underestimated, this dis- 
crepancy is not overwhelmingly large. Interestingly, Doppler 
images commonly also show hot surface features (cf. Korho- 
nen et al. 2007; Lindborg et al. 2011; Hackman et al. 2011). 
These may be artefacts of the Doppler imaging procedure, but 



it is not ruled out that they could arise from the anti-cyclonic 
vortices seen in the present study. 

It is obviously very hard to explain the active longitudes and 
their drift based on the vortex-instabiUty scenario; we believe 
that a large-scale dynamo process is responsible for these ba- 
sic features, as commonly believed (e.g. Krause & Radler 
1980; Moss et al. 1995; Tuominen et al. 2002). Neverthe- 
less, it is possible that the vortex-instability contributes to the 
formation of starspots, and may interfere with the dynamo- 
instability, especially during the epochs of lower magnetic 
activity of the stellar cycle. Although it is very hard to pre- 
dict the implications of the vortices in the magnetohydrody- 
namic regime, it would appear natural that spots, either cool or 
warm, generated by a hydrodynamic vortex-instability, could 
also contribute to the apparent decorrelation of magnetic field 
from the temperature structures. 

The influence of the cyclones and anticyclones on the net 
helicity, important for the amplification of the magnetic field, 
is either close to zero (cyclones) or to decrease the net helic- 
ity (anti-cyclones). This would imply that the magnetic field 
amplification would be equally or even more difficult in the 
regions of the vortices; this picture, however, may be compU- 
cated by the presence of strong horizontal motions present in 
these structures, that might amplify the magnetic field simply 
by their capability for advecting the field lines. 

4. CONCLUSIONS 

We report the formation of large-scale vortices in rapidly 
rotating turbulent convection in local f-plane simulations. The 
vortices appear provided the Reynolds and Coriolis numbers 
exceed critical values. Near the critical Coriolis number, 
the vortices are cyclonic and cool in comparison to the sur- 
rounding atmosphere, whereas for faster rotation warm anti- 
cyclonic vortices appear (see also Chan 2007). The relative 
temperature difference between the vortex and its surround- 
ings is of the order of five per cent in all cases. This is of 
the order of the contrast deduced indirectly from photometric 
and spectroscopic observations of late-type stars. In our sim- 
ulations the typical size of the vortices is comparable to the 
depth of the convectively unstable layer. However, we have 
not studied how the size of the structures depends e.g. on the 
depth of the convection zone. 

We propose that the vortices studied here can be present in 
the atmospheres of rapidly rotating late-type stars, thus con- 
tributing to rotationally modulated variations in the brightness 
and spectrum of the star. Such features have generally been 
interpreted to be caused by magnetic spots, reminiscent to 
sunspots. However, our results suggest that the turbulent con- 
vection and rapid rotation of these stars can generate large- 
scale temperature anomalies in their atmospheres via a purely 
hydrodynamical process. Similar vortex- structures are ob- 
served in the atmospheres of Jupiter and Saturn. Although 
their definitive explanation is still debated, it is possible that 
they are related to rapidly rotating thermal convection in the 
atmosphere. 

However, several issues remain to be sorted out before the 
reality of cyclones and anti-cyclones in the surface layers of 
stars can be established. The current model is highly simpU- 
fied and neglects the effects of sphericity and magnetic fields. 
In spherical geometry more realistic large-scale flows can oc- 
cur which might lead to other hydrodynamical instabilities. 
However, current rapidly rotating simulations in spherical co- 
ordinates have not shown evidence of large-scale vortices (e.g. 
Brown et al. 2008; Kapyla et al. 2010, 2011b), altiiough non- 
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axsisymmetric features are seen near the equator (Brown et al. 
2008). It is possible that the lack of large-scale vortices in 
these simulations is related either to the lack of spatial resolu- 
tion or too short integration time. 

Magnetic fields, on the other hand, are ubiquitous in stars 
with convection zones. Furthermore, on the Sun they form 
strong flux concentrations, i.e. sunspots. At the moment, di- 
rect simulations cannot self-consistently produce sunspot-like 
structures in local geometry (e.g. Kapyla et al. 201 la). How- 
ever, the magnetic fields in global simulations are also very 
different from the high-latitude active longitudes deduced 
from observations, namely showing more axisymmetric fields 
residing also near the equator (e.g. Kapyla et al. 2010; Brown 
et al. 2011). The apparently poor correlation between mag- 
netic fields and temperature anomalies in surface maps based 
on Doppler imaging also suggests that an alternative mech- 
anism might be involved. The presence of large-scale high- 
latitude vortices presents such an alternative. 

Currently it is not clear what happens to the vortices when 
magnetic fields are present. Our previous dynamo simula- 
tions in the same parameter regime (Kapyla et al. 2009) did 
not show clear signs of vortices in the saturated regime of the 
dynamo although this might be explained by the too short in- 
tegration time. Addressing this issue, however, is not within 
the scope of the present paper and we will revisit it in a future 
publication. 

The simulations were performed using the supercomputers 
hosted by CSC - IT Center for Science Ltd. in Espoo, Finland, 
who are administered by the Fiimish Ministry of Education. 
Financial support from the Academy of Finland grants No. 
136189, 140970 (PJK) and 218159, 141017 (MJM), and the 
'Active Suns' research project at University of Helsinki (TH) 
is acknowledged. The authors acknowledge the hospitality of 
NORDITA during their visits. 
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